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ABSTRACT 


In this paper, we examine the spacial distribution of gamma-ray bursts 
(GRBs) using a sample of 373 objects. We subdivide the GRB data into 
two redshift intervals over the redshift range 0 < z < 6.7. We measure the 
two-point correlation function (2PGF), ^(r) of the GRBs. In determining the 
separation distance of the GRB pairs, we consider two representative cosmo¬ 
logical models: a cold dark matter universe plus a cosmological constant A, 
with (r2jn,r2A) = (0.28,0.72) and an Einstein-de Sitter (EdS) universe, with 
(Gm, = (1)0). We hnd a 2 ;-decreasing correlation of the GRB distribution, 
which is in agreement with the predictions of the current structure formation 
theory. We ht a power-law model ^(r) = (r/ro)”'*' to the measured ^(r) and ob¬ 
tain an amplitude and slope of ro = 1235.2 ± 342.6 h~^ Mpc and 7 = 0.80 ± 0.19 
(la conhdence level) over the scales r = 200 to 10^ h~^ Mpc. Our result provide 
a supplement to the measurement of matter correlation on large scales, while the 
matter distribution below 200 h~^ Mpc is usually described by the correlation 
function of galaxies. 


Subject headings: observations - gamma-ray bursts 
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Introduction 


Redshift surveys of galaxies have long been carried out to study the large-scale 
structure (LSS) of our Universe. The two-point correlation function (2PCF) is a statistic 
that can be easily determined from a well-measured galaxy sample. Theoretical predictions 
of the 2PCF can also be made from different dark matter models and structure formation 


scenarios. Thus the 2PCF of galaxies or quasi-stellar objects (QSOs, or q uasars) 


las been 


used as an important statistic to disti nguish different theoretical models flKundic 


Eisenstein et ah 


20051: 


Ross et al 


20071) 


19971: 


Gamma-ray bursts (GRBs) are the most energetic events known to occur in the 
Universe. They are often associated with the death of massive stars and have a redshift 
up to ~ 8. They are robust indicators of the matter-dense part of the intermediate- or 
even the high-redshift Universe. Their spacial distribution helps explore the LSS of the 


Universe. Although recent works suggested t 


in the sky fiGai et al 


2013 


Ghang et al. 


l at GRBs have an ani sotropic distribution 


2014 : 


Wang fc Wane 


demonstrated that GRBs are distributed isot ropically flBriggs 


2014r), most studies 


have 


Magliocchetti et al 


2003 


Vavrek et al. 


19961: 


Meszaros et al. 


2000 


20081) . the latter of which is demanded by the 


cosmological principle. 

Besides its isotropic distribution in the sky, a homogeneous distribution of the GRBs 


is also expected. The d iscovery of a possib 


been recently reported flHorvath et al. 


e structure in the GRB sky distribution has 


20141) . It has a redshift 2 ; ~ 2 (at a distance of 


approximately ten billion light years away) and its size is about 2000 to 3000 Mpc. This 
excess clustering of GRBs has a statistical signihcance of 3cr conhdence level (c.l.) and 
therefore cannot be attributed to the sampling biases at this redshift. In the current 
structure formation theory, structures in today’s Universe like galaxies and galaxy clusters 
etc. all stem from the primordial Gaussian random phase fluctuations of the mass density 
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field. Given the finite time from the end of the cosmic inflation to the present, the evolntion 
of the non-linear strnctnres in onr Universe is limited. They shonld not be larger than 
the scale on which the density contrast beco mes <5^ ^ 1 at present as predicted by the 


linear growth theory of density pertnrbationqj flLabini fc Vasilyev 


2008 1. Since GRBs are 


potential tracers of normal matter, the discovery of this GRB strnctnre (if conhrmed by 
fntnre investigations) casts new shadows on the cosmological principle as well as the cnrrent 
strnctnre formation theory. In fact, a nnmber of interesting resnlt s have been report ed 


on the homogeneity scale of the galaxy and qnasar distribu tions flGroom et al 


2005 


Yadav et al. 


2005 


Sarkar et al. 


2009 


Scrimgeour et al. 


201211 . To do a similar clustering 


analysis on the distribution of GRBs, the correlation function of the GRB distribution has 
to be hrst measured. 


In this work, we measure the real-spac e 


We use the catalogue presented by 


two-p oint correlation function ^(r) of GRBs. 


Greiner! (1201411 . It contains 373 objects over a redshift 


range of 0 < 2 ; < 6.7 by September 3rd, 2014. We subdivide the samples into two redshift 
regimes: z < 2 and z > 2. The 2PGF of each regime is calculated. We ht a power law to 
both of the measured ^{r). In numerical analysis, we prefer the AGDM (cold dark matter 
plus a cosmological constant) cosmological model, with (flni,f^A) = (0.28,0.72) over the 
Einstein-de Sitter (EdS) Universe, with (flm, Ha) = (1)0). 


The rest of the paper is organized as follow: in Section 2, we introduce the methodology 
of the estimation of the correlation function, including the estimators, the biases, and the 
error calculations. In Section 3.1, we describe the data we use and the auxiliary random 
sample. In Section 3.2, we plot the measured ^(r) for the GRBs. In Section 3.3, we ht 
a power law to the measured ^(r) and give the best-ht values of amplitude and slope. 


^At the present time, on a scale smaller than Tc, the growth of density contrast would 


have become nonlinear. 
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Conclusion and discussions are presented in Section 4. 


Estimation of the correlation function 


2.1. The two-point correlation function and the estimators 


The 2PCF, is defined by the probability of finding an object in the spacial volume 


dVi with another object in that is separated by a distance x, i.e. fjPeebles 

dP = n^[l + ^{x)]dVidV2. 


19801 


( 1 ) 


n is the mean number density of the objects. For statistical estimation of ■^( 2^)5 auxiliary 
random sample of Nji points must be generated in a window W. A window IF is a 


three-dim ensional space o 


Following 


Kerscher et ah 


volu me V equivalent to that on which the observation is made. 


20001) . we define the pair count with a unitary function F(x, y): 

( 2 ) 


Nonir) = ^^F(x,y). 

xGD yeR 

The summation runs over all the coordinates of objects (represented by x and y) in the 
observed data set D and the random sample R in the window W. The value of the 
function E(x, y) equals 1 when the separation of the two objects is within the distance 
d(x, y) e [r — Ar/2,r -|- Ar/2] or otherwise equals 0. d(x, y) is the comoving distance 
between the two objects and Ar is the bin width being used in the statistical estimation of 
^{r). With these preparations, we can determine the 2PCF from the observed data. 


Several estimators of ^{r) are popular. One was presented by 
i.e. the DP estimator, 

“ MW “ 


Davis fc Peebles! ( 19831 ). 


( 3 ) 


DD{r) is the normalized number of all pairs of GRBs in the observed data set and is given 
as DD{r) = Niju{r)/[Nj:,[Nij — 1)], where is the total number of GRBs in the data 
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set. DR{r) = NoR{r)/{N^Nji) is the normalized number of pairs between the data and the 
random sample with separation of (i(x, y) G [r — Ar/2,r + Ar/2]. NddIt) is dehned in a 
similar way as Nuji{r) in equation ([2]), while ‘D’ and ‘i?’ refer respectively to the observed 


data set of GRBs and the auxiliary random catalogue. Two o 


are the Landy-Sza’ 
estimator ^Ham(?^) I 


ay fLS) estima tor .^ls(’") flLandy fc Szalav 


her w idely used estimators 


Hamilton 


1993) and the Hamilton 


199311 : 


4s (r) = 

C^Ham ) 


DD{r) — 2DR{r) + RR{r) 
RR{r) 

DD{r)RR{r) ^ 

" [DR{r)Y 


( 4 ) 

( 5 ) 


RR{r) = N/[Nji{NR — 1)] refers to the normalized number of pairs with the separation 
mentioned above in the random sample. is dehned in a similar way as N£,fi{r) in 

equation ([2]). 


2.2. The biases 


A relevant problem for estimating the 2PCF is that there might be systematic biases 
and stochastic noise which perturbs any real determination of the 2PCF, especially on 
small-amplitude values of At small distances, the estimators mentioned above have 
very similar performances. However, on large scales, they are not totally equivalent and 
some of them could be biased. 


Generally speaking, for a valid statistical estimator of the quantity X, the sample 


average A in a hnite volume V, i.e. X{V), must satisfy (IGabrielli et al.l 120041) 


lim X{V) 

V—>-oo 


w, 


( 6 ) 


where (X) is the ensemble average. For a hnite sample volume V, there is a systematic bias 


in the estimator X. An unbiased estimator is the one that satishes X{V) = (X). 
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The LS estimator is generally biased as the DP and Hamilton estimators. There are 
also other two biased estimators that are commonly used, i.e. the natural estimator and 
the Hewett estimator. Nu merical tests on artihcial distributions have been d one to study 


biases in these estimators (jPons-Borderia et ah 


1999 


Labini fc Vasilyev 


20081). It has been 


found that on large scales, the LS estimator a nd the Hamilton estimator signihcantly 


outperform the rest (jPons-Borderia et ah 


19991) . The LS estimator has an indistinguishable 


performance as the Hamilton estimator but is less sensitive to the size of the random 
sample, Nj^. Thus it is more preferable from a practical point of view. 

In this work, we focus on the LS estimator of while the DP and Hamilton 
estimator ar e also used for c ompa rison. A detailed comparison of these estimators can 


be found in 


Kerscher et ah 


fl2000[l . while alternative estimator s (such as the f ull-sh ell 


estimator and the geome tric estimator, etc.) were proposed in 
Labini fc Vasilyev ( 2008 1. 


Kerscher et ah 


(1200011 and 


2.3. The errors 


A limited sample would result in a large likely error or variance. A pr^er way of 
estimating errors could reduce this effect to the least. There are several wayo to determine 
the measured errors of the 2PCF of GRBs. Three most common ones are the Poisson 
estimate, the ‘£eld-to-£eld’ (FtF) error, and the jackknife estimate. The Poisson estimate 
of the errors of ^(r) is given as: 

1 + ^(r) 


o-poi(r’) = 


DD(r) 


(7) 


^One can refer to 


has also been given in 


Hamilton! (1199311 and see the r eferences therein. A so-called ‘cox process’ 


Pons-Borderia et ah 


(1199911 . 











































The second method is the ‘FtF’ method. In this method, the whole sample is 
divided into Ni, snbsamples. The 2PCF of each snbsample is calcnlated, i.e. ^i{r) with 


i = 1, 2,iVft. The error is calcnlated by 

Nb 


a- 


FtF 


r = 


1 V—^ DRAF . n 


AT, - 1 ^ DR{r) 


( 8 ) 


^{r) is the estimate of the 2PCF on the entire sample. For onr stndies, the entire sample of 
GRBs, which has a redshift range of z G [0, 6.7], is divided into seven redshift bins with the 
interval of 1. Thns there is = 7. Details of the snbsamples are given in Table [B 

The third way to estimate the errors of ^(r) is called the jackknife estimate. It is given 


as 


(Ti 


r = 


N' 

E 


DRiir) 




(9) 


^ DR{r) 

Like the ‘FtF’ method, the entire sample is divided into N' snbsamples. denotes the 

estimate of the 2PCF on all of the {N' — 1) snbsamples except the i-th one. 

On small scales r < Mpc {h is the Hnbble constant in nnits of 100 km s“^ 

Mpc“^), all these errors have comparable magnitndes, while above this scale (np to 
r ~ 100h ~^ Mpc), the jac kknife and ‘FtF’ errors are considerably larger than the Poisson 


estimate flRoss et al 


2007(1 . In addition, on small s cales, the jackknife method gives smalles t 


flnctnations of the three and is thns more efficient fjScranton et al. 


2002 


Zehavi et al. 


2OO2I). 


So in this work, we make nse of the jackknife method to estimate the errors of the measnred 
2PCF. Like the ‘FtF’ error, we divide the whole GRB sample into 7 snbsamples as well as 
7 redshift bins. Each snbsample corresponds to one of these redshift bins. Thns there is 
N' = 7. 


















3. Data and Results 


3.1. The GRB data and the random catalogue 


We use the GRB catalogue! presented by J. Greiner to determine the 2PGF. It is a 
subjective collection of GRBs that are detected by a number of satellites and programs, 
i.e. BATSE, RXTE, BeppoSAX, HETE-2, Interplanetary Gamma-Ray Burst Timing 
Network (IPN), INTEGRAL, Swift, AGILE, and the Fermi Gamma-Ray Observatory. The 
catalogue contains thousands of objects and is updated almost everyday. We use the data 
that were released before September 3rd in 2014. Among those 398 objects which have the 
redshift measured, only 375 of them are well determined. 300 (about 80%) of them come 
from one source, the Swift satellite. The rest of the samples come from Fermi Gamma-Ray 
Observatory, HETE-2, and other telescopes. For the others only an upper limit is provided. 
We subdivide these 375 GRBs into 10 redshift intervals with effective redshifts from z = 0 
to 9.2. In the redshift bins with z > 7 in the catalogue, there are only two GRBs: the GRB 
090423 {z = 8.26) and the GRB 090429B {z = 9.2). They are omitted in our studies since 
they have little statistical signihcance. 


Therefore, we base our research on these = 373 GRBs that have redshifts and 
angular positions well determined. They cover a redshift range of 0 < 2 ; < 6.7 and we use 
N' = 7 for the estimates of ^(r) and the jackknife errors. Robustness of the estimators 


introduced in t 

le last section has been well c 

emonstrated for a sample of hundreds of 

obiects ( 

Shaver 

1964, 

1967; 

Shanks et al. 

1987 

Anderson et al. 

1968; 

lovino & Shaver 

1988 

Mo & Fane 

199 

3 


Shanks & Bovle 

1994 


Groom & ShankJ 

1996 

: Pons-Borderia et al. 

I999I; 


Kerscher et al 


2000]). The GRB data we use are listed in Table 3 which is publicly available 


^ See http://www.mpe.mpg.de/~jcg/grbgen.html 
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onlinqj. The redshift distribution of the data are presented in Tableland Figure [T] The 
density of random points used is 20 times the density of GRBs. We use bin widths of 
Alog(s) = 2h~^ Mpc. 


3.2. The real-space two-point correlation function ^(r) 


The measured values of 2PCF heavily rely on the assumed cosmology in that the 
comoving separation of the GRBs are model dependent. In this work, two cosmological 
models are considered: 1) the Lambda-cold dark matter (AGDM) cosmological model which 


uses the parameters det ermined from 


he nine-year Wilkinson ]V 


2013 


Hinshaw et ah 


20i; 


icrowave Anisotropy Probe 
) in conjunction with other 


(WMAP) observations fiBennett et ah 

observation j^, with (Gm, Ga) = (0.28,0.72); 2) an Einstein-de Sitter (EdS) cosmological 
model, with (r2jn,r2A) = (1)0). 

The 2PGFs of the GRB data for different densities of the random points are presented 
in Figure [21 We plot these to study the potential impacts of the statistical noise on the 
results. They cover a distance scale up to ~ 800h“^ Mpc. Most GRB pairs in the data have 
a separation distance over 100h“^ Mpc. The current structure formation theory predicts 


fISpringel et al. 


2005 


the clustering o 

' matter remains in 

Eisenstein et al. 

2007 

). The redshift 


The GRBs data we use is presented in a text format, which is publicly available online. 


® Other observations include 


spectrum data (IFowler et al. 


2010 : 


:he high-£ cosmic microwave bac 


Das et al, 


the baryonic acoustic oscillation (BAO 


Anderson et al. 


2012 


Blake et a. 


® See Section 9.4 in 


Dodelsonl (120081) 


aground (GMB) 


powe r 


2012), and a new Hq measurement (IRiess et al 


2011 

Keisler e 

al. 

on 

Reichardt et al. 

2012) 

(Bentler et al. 

2011 

Padmanabhan et al. 

2012 


201111 . 
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small-scale pe culiar velocities of the objects and the redshift variances are also minimal 


on this scale flRoss et al 


20071) . Thus the difference between the redshift-space and the 


real-space correlation functions on such large scales could be negligible. 


Two other comments on Figure [2] are necessary. The hrst one is that the ^(r) of GRBs 
is quite scattered on scales below 100h“^ Mpc. The poor performance of the estimates of 
,^(r) and its errors on such a scale is due to the lack of data for r < 100h“^ Mpc. Most 
GRBs have been found to have a separation with each other above 100h“^ Mpc. More 
observations are needed to improve the performance of the analysis on r < 100h“^ Mpc. On 
the scale r > 100h“^ Mpc, the results don’t show much difference for varying the objects 
density in the random catalogue. Thus they are more reliable and can be used to hnd the 
possible power law ht of ^(r). 


The second one is that it shows a ‘bump’ at r = 200h.“^ Mpc in the measured .^(r) 
for the GRB sample. For an EdS Universe, it is located at r ~ 100h“^ Mpc. For different 
number densities of the random points, the bump exists, implying that it cannot be 
attributed to statistical noise. In cosmology, the temperature drop of the Universe at about 


380,000 years after the big bang would cause a sudden declin e of sound speec 


fluid. This would leave the oscillations in it become frozen (IMartinez et al. 


in the matter 


2009). For 


the AGDM model with (GuijGa) = (0.28,0.72), the signal manifest s as a peak at about 
100h“^ Mpc in the correlation function for the galaxy distribution flEisenstein et al.ll2005l) . 


which is about 100h“^ Mpc away from the ‘bump’ we discovered. However, since the size 
of the GRB sample we use is not sufficiently large, whether these two have any physical 
connections or not is still subject to future investigations. 


To study the redshift evolution of the 2PGF for GRBs, the entire GRB samples are 
subdived into two groups with z < 1.5 (173 objects) and 2 ; > 1.5 (200 objects). The 2PGFs 
are calculated for each group for the AGDM model with (Dm, Da) = (0.28, 0.72) and are 

















12 


plotted in Figure |3]for comparison. From Figure |3l it is found that on the scales r < 

Mpc, the low- 2 ; group of GRB samples have a higher correlation amplitude than the high- 2 ; 
group. This is compatible with the current structure formation theory predicting that the 
precursors of GRBs with low redshifts would have more time to grow and thus would become 
more correlated with each other. The differences become small and indistinguishable on 
scales r > 5 X Mpc. In fact, some previous stud i es have suggested a z - increasing 


corre la tion behavior in quasars and halos ([Kundk 


200511 . 


1997 : 


La Franca et al. 


1998 : 


Gao et al. 


Knndid fjl997l) had reported that the amplitude of the quasar correlation function 


of the high-redshift samples {z > 2 ) was significantly higher than that of the low-redshift 
sample {z < 2). As shown in Figure |3l this is not the case for the GRB samples. More 
discussions are presented in Section 01 


3.3. A power law fit to ^(r) 


A power law of the form 


e(r) = - 

To 


-7 


( 10 ) 


has been respectively htted to the 2PGFs of galaxies and galaxy clusters for a cold dark 
matter (CDM) Universe (flm = 1), using a least-y^ technique. For the case of galaxies, on 
scales r < 10/i“^ Mpc, the correlation length rn has a best-ht yalue of 3.76 < rn < 7.3 h ~^ 


Mdc, with th 

e exponent index 7 

.5 < 7 

Hermit 

1996; 

Zehavi et al. 

2002 , 

2004) 


. 5 < 7 < 1.8 flDavis fc Peebles 


1983 : 


Maddox et al 


1990; 


20041) . The values vary somewhat from literature to 


literature depending on the galaxy sample utilized, the estimator used, and the weighting 


scheme employed, etc. For the case of galaxy cluste rs, the 2 


i.e ^rr(r) = fr/rn) with 13 < 

rn < 25 h-^ 

Vlpc 

Po.stman. Hnchra. & Geller 

1992; 

Nichol et al. 

1992 


'CF a 


1988 


Dalton et al. 


so follows a power law . 


Peacock fc West 


1992 . 


199J). 


1992 
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A power-law model of the same form as equation fllOl) has also been htted to the 
measured .^(r) of GRBs above 100/i“^ Mpc in Figure 01 The parameters tq and 7 are 
determined by a least-y^ procedure. In the ACDM cosmology, we hud a best-£t power law 
with To = (1235.2 ± 342.6) h~^ Mpc and 7 = 0.80 ± 0.19 (la conhdence level) on scales 
r = 200 to lO'^ h~^ Mpc. In the EdS Universe, the results are ro = (322.4 ± 92.3) h~^ Mpc 
and 7 = 0.62 ± 0.20 (la conhdence level). The values of tq and 7 for both cosmological 
models are presented in Table [2l The scale-length ro in the power law ,^(r) = (r/ro)”"^ is 
not very well determined in the analysis. That is because the data points of measured ^ 
have large error bars, due to the limited size of the GRB sample. This would be improved 
once a larger sample of GRBs with well-determined redshifts are available for the analysis. 

Most galaxies have a redshift z < 1.5. For a comparison with the best-£t power law to 
that of galaxies, we £t a power law to those GRBs with z < 1.5. This sub-catalogue contains 
173 objects. We obtained a best-fit rg = 181.4 ± 113.9 h~^ Mpc and 7 = 1.12 ± 0.46, with 
Xmin = 0.54. The best-£t result was presented in Figure [5l The la-error bars are calculated 
from the jackknife method by equation ([9]) with the N' = 5. 


Conclusions and Discussions 


Up to now, the studies of the correlation function ^{r) have usually been limited to 
the galaxy samples at low redshift, i.e. 0 < 2 ; < 1. In this paper, we extended this work 
to the GRBs samples. Many of them have a redshift range oi z > 1 and thus can be used 
to explore the matter clustering and its ev olution in the earlier Universe. We calculated 


the 2PGF ^{r) for the GRB samples from 


Greiner! (120141 ). We considered two popular 


cosmological model in estimating ^(r), a AGDM cosmology with (flm^f^A) = (0.28,0.72), 
and an EdS one with (flnn ^^a) = (1)0). We found that in the AGDM Universe, the .^(r) 
for the GRB samples on large r (~ 200 to 10^ h~^ Mpc) follows a best-fit power law with 
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ro = (1235.2 ± 342.6) h~^ Mpc and 7 = 0.80 ± 0.19 (la confidence level). For the EdS 
Universe, the results are tq = (322.4 ± 92.3) h~^ Mpc and 7 = 0.62 ± 0.20 (la conhdence 
level). Both correlation functions showed systematic deviations from the best-£t power 
law below the scale r = 200h“^ Mpc. We concluded these departures are due to the 
volume-limited sample, which has insufficient GRB pairs with separation below 100h“^ 
Mpc. In fact, on the scale r < 100h“^ Mpc, the matter distribution in our Universe is 
usually determined by the redshift galaxy survey, such as the Sloan Digital Sky Survey 


(SDSS) f Aihara et al. 


201111 . Our work provided a supplement for the measurement of the 
matter correlation on scales r > 200h“^ Mpc. For a better performance of the measurement 
of ^(r) for GRBs on smaller scales, which can be used to compare with the results obtained 
from the galaxy survey below 100h“^ Mpc, further observations are needed. 

In fact, since the GRB sample we used is not large enough, one has to take 


these resul 


s with a grain of s alt. Besides, the Wilkinson Microwave Aniso tropy Probe 


(WMAP) ( Bennett et al 


2013 1 and the PLANGK satellite ( Ade et al. 2013 1 have provided 


unprecedentedly precise measurements of the anisotropy and the inhomogeneity of the 
matter distribution in the early Universe. To answer the questions mentioned above, one 
may have to use the WMAP/PLANGK observational data for a combined analysis. A more 
careful clustering analysis of the medium redshift Universe using the WMAP/PLANGK 
observations together with the quasars and GRB data is currently undertaking. We 
hope that the results would shed new light on the structure formation theory and the 
inhomogeneities of our early Universe. 


We are grateful to Zhi-Bing Li from the Sun Yat-Sen University and Zhe Ghang 
from the Institute of High Energy Physics for all the useful discussions and suggestions 
The final version of the manuscript has greatly benefited from the careful scrutiny of 
the referee. This work is based on the GRB catalogue presented by Jochen Greiner at 
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Fig. 1.— The redshift distribution of the 373 GRB data. The y-axis denotes the number of 
objects in each redshift bin. 
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Fig. 2.— The real-space correlation function, ,^(r) of the GRBs for different number of 
objects in the auxiliary random catalogue. The results are presented in both the ACDM 
model, (r2m,r2A) = (0.28,0.72) and the EdS Universe, (r2m,r2A) = (1,0). Circles, squares, 
and triangles respectively represent the results obtained from that the random catalogues 
are 5, 10, and 20 times the density of GRB data. Error bars are calculated from the jackknife 
method by equation ([9]) with the N' = 7. 
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Fig. 3.— The redshift evolution of the measured .^(r) for GRBs in the ACDM model with 
= (0.28,0.72). The ^{r) for 2 ; < 1.5 (black circles) and for 2 ; > 1.5 (blue squares) 
with IcT jackknife error bars are plotted. 
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(a) (f2„„f2A) = (0.28,0.72) (b) = (1,0) 

Fig. 4.— The best-fit power law of the measured .^(r) for the GRBs from iGreinerl fl2014j) . (a 


For the AGDM universe with (him, ^a) = (0.28,0.72). The best-£t values of the correlation 
length ro and the slope 7 are ro = 1235.2 ± 342.6 h~^ Mpc and 7 = 0.80 ± 0.19 (la 
conhdence level), with = 0.19. (b) For the EdS Universe with (him, Ha) = (1,0). The 
best-£t parameters are ro = 322.4 ±92.3 h~^ Mpc and 7 = 0.62 ±0.20 (la conhdence level), 
with Xmin ~ 0.21. In both hgures, triangles represent the measured .^(r) obtained from that 
the density of random points is 20 times of that of the GRB data. The la-error bars are 
calculated from the jackknife method by equation ([9]) with the N' = 7. 

























24 



Fig. 5. - 


Th e best-fit power law of the measured .^(r) for the GRBs at z < 1.5 from 


Greinerl 1 20141) for the AGDM universe with (flm, Q\) = (0.28, 0.72). Triangles represent the 


measured ^(r) obtained from that the density of random points is 20 times of that of the 
GRB data. The la-error bars are calculated from the jackknife method by equation (|9]) with 
the N' = 5. Solid black lines indicate the best-ht power law of ^{r). The best-£t values of 
the correlation length and the slope 7 are tq = 181.4±113.9 h ^ Mpc and 7 = 1.12T0.46, 
with Xmin = 0-54. 
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Table 1: Properties of the GRB catalogue presented by 


Greinerl (120141 ): the redshift bins and 


the number of GRBs in each of the Nh = 7 bins that are used to calculate the ‘FtF’ error. 


Redshift 

Number 

0< z<l 

112 

1 < z<2 

108 

2< z<3 

80 

3 < z < 4 

40 

4 < < 5 

21 

5 < ^ < 6 

8 

6 < z < 7 

4 


able 2: The power-law £t to the measured .^(r) of the GRBs catalogue presented by 


Greiner 


(120141) . The best-£t values of the correlation length tq and the slope 7 for a ACDM universe 


with (flm, = (0.28,0.72) and an EdS universe with = (1)0) over the scales 

0 < r < 800 h~^ Mpc are presented. Xmm = ~ d.o.f + 1) represents the reduced 

minimal chi-squares of the results on the scale r > 100h“^ Mpc. = 9 is the number of the 
bins, d.o.f = 2 represents the degree of freedom, i.e. the number of free parameters in the 
£t. 



ACDM 

EdS 

ro 

1235.2 ± 342.6 

322.4 ±92.3 

7 

0.80 ±0.19 

0.62 ± 0.20 

Xmin 

0.19 

0.21 


















